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oxide thin films near the superconductor-insulator transition 
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Current-voltage characteristics in the insulator bordering superconductivity in disordered thin 
films exhibit current jumps of several orders of magnitude due to the development of a thermally 
bistable electronic state at very low temperatures. In this high-resolution study we find that the 
jumps can be composed of many (up to 100) smaller jumps that appear to be random. This indicates 
that inhomogeneity develops near the transition to the insulator and that the current breakdown 
proceed via percolative paths spanning from one electrode to the other. 

PACS numbers: 74.78.-w,72.20.Ht,73.50.Fq,74.81.Bd 



Disordered, thin-film superconductors can be driven 
into an insulating state by varying an experimental pa- 
rameter such as the level of disorder or the magnetic field 
(B) (for a review, see Ref. [T]). The B-driven insulat- 
ing state is nonmonotonic, the resistance showing a pro- 
nounced peak as a function of B [2H5]. At B's much 
beyond the peak the insulating behavior is replaced by a 
weakly insulating, or perhaps metallic [Hl^, state. Sev- 
eral theoretical groups have addressed this nonmonotonic 
insulator [Tlfl^. but a consensual framework is yet to 
emerge. 

One key observation in both the disorder-driven and B- 
driven insulators is that the current- voltage (I-V) charac- 
teristics exhibit a voltage threshold Vth, separating high- 
resistance (HR) and low-resistance (LR) states. At low 
temperatures (T < 0.2 K), Vth is accompanied by a large 
(several orders of magnitude) current jump, as well as 
hysteresis [131 El- This switching behavior was inter- 
preted in light of a model P^HTTj in which Joule heating 
by the current flow can induce a nonequilibrium state 
where the electrons' T, Tei, can signiflcantly exceed the 
phonon bath T. One consequence of this model is that, 
at very low T, the electrons enter a bistable region where, 
as a function of voltage, T^i can abruptly change causing 
the large current jumps between the "cold" HR state and 
the "hot" LR one. However, there have been some indi- 
cations to suggest that this is not the whole picture. In 
recent work on insulating TiN thin films, multiple con- 
ductance steps in the dl/dV curves were observed near 
the threshold at low B's [TH]. This preliminary finding 
was taken as an indicative of a percolative behavior, re- 
flected by multiple charge depinning transitions |19j , with 
the threshold behavior later attributed to heating insta- 
bility [n]. 

In this Rapid Communication we report the results of 
high-resolution I-V measurements in the vicinity of Vth- 
These more detailed studies reveal the existence of several 
(as many as 100) discontinuous / steps, accompanied by 
hysteresis, near Vth- A thorough experimental analysis of 
these multisteps supports a scenario in which the current 
distribution develops significant inhomogeneity near the 



superconductor-insulator transition (SIT): Each / step 
corresponds to a sudden appearance or disappearance of 
a narrow conductive path of "hot" electrons that carries 
relatively high current through the "cold" insulator. 

The system we study consists of highly disordered 
amorphous indium oxide (a:InO) that was fabricated by 
e-gun evaporation of 99.99% pure sintered In203 pieces 
onto oxidized silicon wafers [501 HI] • Premade Au leads 
and a lift-off mask were used to pattern devices with a 
well-defined two-terminal geometry. 25—30 nm thick and 
500—1000 /im long square films were used throughout 
this study. In addition, one hall-bar sample deposited 
through a shadow mask and contacted by hand-pressed 
indium was used to test the influence of the contact re- 
sistances on the results. Some of the samples were heat 
treated at 40 °C in vacuum for a few hours postdeposi- 
tion to reduce their degree of disorder by annealing [21] . 
All measurements were performed using a standard two- 
probe technique inside a dilution refrigerator capable of 
reaching 0.02 K. The I-V curves were traced by sweeping 
dc voltage and monitoring the current with a low-noise 
transimpedance amplifier. The resistance of the devices 
was measured by a low-frequency (0.3—3 Hz) lock-in con- 
figuration, or by the low-voltage derivative of the I-V 
characteristics. 

A typical switching behavior in our samples can be 
seen in Fig. [ija), where we present a T = 70 mK I-V 
curve taken from sample LWla. At low V, I through 
the film is small (< 1 pA), inapparent on the scale of the 
plot. As V reaches a certain voltage, referred to as Vhl, 
conduction starts with / increasing abruptly by several 
orders of magnitude. On reversing the sweep direction, 
hysteresis is observed and the LR state persists until a 
second voltage, Vlh, is reached. The switching behavior 
is even more pronounced when plotted using a logarith- 
mic ordinate [Fig. [l]jb)], better illustrating the range of 
variation in / at the two thresholds. 

The central result of this work is shown in Fig. [ijc). 
Focusing on the data obtained from the negative sweep 
direction close to Vlh reveals a series of discontinuous 
steps. These steps seem random in location, as well as 
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Figure 1. (a) I-V characteristic of sample LWla at T = 70 
mK and B = 2 T. The arrows show the voltage sweep direc- 
tion; the positive sweep is shown in red, while the negative 
sweep is in blue. This color scheme will be followed in the sub- 
sequent figures, (b) Same as (a) plotted on a semilogarithmic 
scale. The threshold voltages Vhl and Vlh are marked by 
the two arrows, (c) Closeup view on the gray rectangle in 
(a). Inset: I-V characteristics of LW3f at T = 30 mK and 
B = 9 T. Each green arrow denotes a turning point of the 
sweep direction. 

jump magnitude, but are highly reproducible after suc- 
cessive scans, even when the two measuring contacts are 
interchanged. The step at the lowest V corresponds to 
the threshold we designated as Vlh, beyond which the 
system reverts to the HR state. 

The / steps are accompanied by a complex hysteresis 
phenomenon. This is manifested in the inset of Fig. [l|c), 
which shows a set of I-V curves obtained by cycling V 
between a voltage higher than the threshold Vhl and a 
voltage just below each step. Several distinct hysteresis 
loops are seen in these data, each associated with a spe- 
cific current step. While this behavior is observed under 
all experimental parameters for the / steps to appear, 
the exact dependence of the subhysteresis on these pa- 
rameters remain to be investigated further. 

Figure [2] shows the T evolution of the / steps measured 
across film LW3f. As T increases, the steps pattern in the 
negative sweep branch deforms gradually, the location of 
each step shifts to lower voltages, and the step magnitude 



Figure 2. I-V isotherms of LWSf at the vicinity of the thresh- 
old, taken at B = 4 T. For clarity, the curves are offset from 
one another by 25 nA. Inset: Expanded V range view of the 
30, 40, and 75 mK scans, showing the shifting of Vhl- 



decreases until the discontinuity vanishes altogether. 

Interestingly, multiple / steps appear also in the posi- 
tive sweep direction (red traces in Fig. |2|. While at very 
low T we observe a single step at Vhl (see the inset of 
Fig. 2 for a wide V range), Vhl rapidly decreases with T 
and, above ^ 50 mK, breaks into multiple steps similar 
to those in the negative sweep direction. Eventually, the 
hysteresis collapses and the two branches merge into one. 

The evolution of the steps pattern as function of B is 
shown in Fig. [3] where we plot, for sample LW3f, the 
dependence of the I-V characteristics on perpendicular 
B. The critical point of the i3-driven SIT for this sample 
is Be = 0.9 T and the magnetoresistance (MR) maxima 
was found to occur at 5 T field. The step profile appears 
only for fields in the range Be < i? < 10 T, beyond which 
the curve is smooth. The location of each step moves to 
higher V's as B is increased, and reaches a peak near 
B = 9 T. Simultaneously, as the steps shift, their height 
increases, reaching a maximum value, and then decreases 
to zero. The corresponding B dependence of Vhl and 
Vlh is given in the inset. As the difference Vhl — Vlh 
approaches zero, near Be and near 10 T, the hysteresis 
closes and multiple steps are observed also in the positive 



3 



700 



600 :^ 



500 - 



400 



300 



200 



100 





116.5 



mean=0.178mV 
std=0.178mV 




117.5 118 118.5 
Voltage (mV) 

25 

20 

>, 
o 

15 




0.5 1 

Steps Separation (mV) 



0.2 0.4 
Step Height (nA) 



64 65 66 

Voltage (mV) 



Figure 4. 1-V curve of LW3d at T = 20 mK and at B = 0. 
The step locations are marked by the gray vertical lines. Inset: 
Wide V range view of the same curve, (b) Histogram show- 
ing the distribution of the V differences between consecutive 
steps. The black dots are of the exponential distribution of 
the same mean, (c) Histogram of the step height distribution. 



Figure 3. (a) I-V curves of LWSf as a function of B at the 
vicinity of the threshold, taken at T = 30 mK. Only negative 
sweep data are shown. For clarity, an offset current of multi- 
ples of 20 iiA has been added. Inset: Vlh and Vhl vs B of 
the same sample. 

sweep direction (data not shown). 

The magnetic field is seen to have only a moderate in- 
fluence on the hierarchy of the step pattern. For each / 
step, a trajectory on the B-V plane can be clearly and 
independently traced from just above Be to its eventual 
high-B demise. Even though these trajectories differ in 
shape, they seldom appeared to cross one another, leav- 
ing the step arrangement almost unchanged over a wide 
range of i?'s. 

We now turn to films with a higher degree of disorder. 
Sample LW3b has a disorder value that is larger than 
critical and is insulating at i? = 0. Nonetheless, it still 
exhibits a MR peak at B = 2 T [T4 1 I22H24 ] . In Fig. ga) 
the 1-V step pattern of this device at T = 20 mK and 
S = is shown. A few dozens / steps can be seen within 
the sensitivity of the measurement, almost an order of 
magnitude more than for the superconductivity samples. 
The step locations, found using an edge detection algo- 
rithm, are delineated by the gray vertical lines. In total, 
92 / steps were detected. 

The large number of steps enables us to carry a statis- 
tical analysis of the process. Using the data in Fig. |4][a). 



we plot two histograms showing the V differences be- 
tween consecutive steps and the step heights. The V dif- 
ference histogram is exponentially distributed, suggest- 
ing that the steps are randomly spread in a Poisson-like 
fashion. The step-height histogram is bell shaped, cen- 
tered on a mean height of 0.17 nA with a low standard 
deviation value of 0.057 nA. For comparison, this anal- 
ysis was performed on the same sample after thermal 
annealing for 72 h (at its postannealing state, this sam- 
ple was referred to above as LWSf). The corresponding 
total number of steps and their mean height, extracted 
from an I-V curve at T = 30 mK and S = 4 T (Fig. [2|, 
is 22 and 1.21 nA, respectively. 

To account for our findings we adopt the perspective 
that transport in these films can become spatially in- 
homogeneous. In this context our results appear quite 
reasonable, considering that the conduction transition is 
not singular and involves multiple events of the same sort, 
randomly arranged, and strongly contingent on the level 
of disorder. Further evidence for such inhomogeneity 
near the SIT have been reported in the past [U [TH [251 - 
and has even provided a basis for several models 
describing the phase transition [OHTTl [T9l [30l [31] . If the 
fluctuations in conductivity are strong, the current flows 
along narrow, percolating, channels of least resistance. 
Combining this view with the electron heating scenario 
from Refs. [111 HHI gives rise to the possibility of local 
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overheating of electrons at each path. 

Within this physical picture, we can attribute the / 
steps to spatially separated bistabilities in T^i. Depend- 
ing on the applied V, the conduction regime of each chan- 
nel can be toggled almost individually between the "cold" 
HR state and the "hot" LR one. The low-T behavior 
can be understood as follows: At low V, all channels are 
nearly closed. When Vhl is reached, conduction begins 
with the opening of all channels in an avalanche fash- 
ion; the first channel to open triggers the conduction in 
nearby channels, which then proliferates to the bulk of 
the sample. In the I-V curve, this process appears as one 
big current step. Upon reducing V, channels close one 
by one, and so only small steps are found in the oppo- 
site sweep direction. This behavior varies somewhat at 
high T, where, as depicted in Fig. [2] we find that small 
steps also appear in a positive sweep direction, indicat- 
ing that the high-T rupture of the "hot" phase is also 
successional. 

From this interpretation, it follows that the / step pat- 
tern serves as a fingerprint of the particular configuration 
of inhomogeneitics. Varying the degree of disorder con- 
siderably influences the low-T inhomogeneity: At lower 
disorder fewer channels are available, but a higher current 
flows through each one. On the other hand, the hierarchy 
of the steps, being largely conserved over a wide range 
of B's, suggests that inhomogeneity is only weakly de- 
pendent on B. This is intriguing, considering the strong 
fleld dependence of the threshold, and certainly calls for 
a more thorough theoretical analysis. 

It is interesting to note the similarities between insu- 
lating a:InO and polycrystalline TiN bordering the SIT. 
Shared key features include the Arrhenius-activated be- 
havior, the pronounced peak in the MR [3l |4j, and the 
high-fleld saturation [5j E] . When comparing the present 
data with those from Ref. [TH], additional distinct sim- 
ilarities are observed. Both materials exhibit a set of 
sharp steplike V thresholds that are nonmonotonic in B 
and attain their maximal value far above the MR peak, 
further suggesting a common underlying mechanism. 

In summary, the appearance of multiple / steps near 
the threshold for conduction can be linked to an inho- 
mogeneous conduction state that develops in the vicinity 
of the SIT. We note that features similar in appearance 
to those reported here have been observed in the context 
of two-dimensional quantum dot arrays |32j . a system 
in which current is known to be inhomogeneously dis- 
tributed due to disorder [33]. These similarities further 
advocate the relevance of Coulomb blockade in the SIT 
insulating phase US] , and prompt the use of quantum 
dot arrays as a model system for understanding its trans- 
port properties. 

We wish to thank V. Kravtsov, Y. Meir, and B. Sacepe 
for illuminating discussions. This work was supported by 
the Minerva foundation with funding from the Federal 
German Ministry for Education and Research. 
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